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Abstract 

' Criteria are formulated both for the existence and for the non-existence of 
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complex eigenvalues for a class of non self-adjoint operators in Hilbert space 
invarariant under a particular discrete symmetry. Applications to the PT- 
symmetric Schrodinger operators are discussed. 



1 Introduction and statement of the results 

(N 

in 

The Schrodinger operators invariant under the combined application of a reflection 



symmetry operator P and of the (antilinear) complex conjugation operation T are 
called PT-symmetric. A standard class of such operators has the form H = H + iW 
Q_i' where: 

1. Hq is a self-adjoint realization of — A+V on some Hilbert space L 2 (f2); f2 C 1Z n , 
n > 1; V and W are real multiplication operators. 

2. V are even and odd with respect to P, respectively: PV = V, PW = —W. P 



is the parity operation 

(Plp)(x) = . . . , (-iy"X n ), 4> G L 2 

where ji = 0, 1; % = 1 for at least one 1 < i < n; 

If T is the involution defined by complex conjugation: (Tip)(x) = ip(x), one imme- 
diately checks that (PT)H = H(PT). 
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PT— symmetric quantum mechanics (see e.g. [I] , [2] , [3] , [I] , |H] |Hj |lj jHj,[?]) requires 
the reality of the spectrum of PT— symmetric operators, recently proved, for in- 
stance, for the one dimensional odd anharmonic oscillators [13], ^2]- Imposing 
boundary conditions along complex directions, however, examples of PT— symme- 
tric operators with complex eigenvalues have been constructed |14j . It is therefore 
an important issue in this context to determine whether or not the spectrum of PT- 
symmetric Schrodinger operators with standard L? boundary conditions at infinity 
is real. We deal with this problem only in perturbation theory, but we will obtain 
criteria both for existence of complex eigenvalues (Theorem 1.1) and for the reality 
of the spectrum (Theorem 1.2), in even greater generality than the PT symmetry. 

Let Ti be a Hilbert space with scalar product denoted (x\y), and H : H — > H 
be a closed operator with dense domain T> C 7i. Let H\ be an operator in Ti with 
T>(Hi) D T>. This entails that H\ is bounded relative to Hq, i.e. there exist b > 0, 
a > such that H-ffi^U < frHifo^ll + G V. We can therefore define on V 

the operator family H e := H t = H + eHi, Ve G C. 

We assume the following symmetry properties: there exists a unitary involution 
J : H — > Ti mapping V to V, such that 

JH = H*J, JH l = H{J (1.1) 

In other words, J intertwines Hq and H\ with the corresponding adjoint operators. 
Note that: 

1. The properties J 2 = 1 (involution) and J* = J" 1 (unitarity) entail J* = J, 
i.e. self-adjointness of J; 

2. The properties (jl.l|) entail, if e G 1Z, JH e = H*J; therefore the spectrum 
cr(H t ) of H e is symmetric with respect to the real axis if e G 1Z. 

3. An example of J is the parity operator P. 

Let Hq admit a real isolated eigenvalue Ao of multiplicity 2 (both algebraic and geo- 
metric, i.e. we assume absence of Jordan blocks). Let e\, be linearly independent 
eigenvectors, and E\ Q the eigenspace spanned by ei,e 2 . Clearly JE Xo '■— E^ o is the 
eigenspace of Hq corresponding to the eigenvalue A = A , and hence the bilinear 
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form (u*\v),u* G El o ,v G E\ is non degenerate. Therefore we can choose ei,e2 in 
E\ in such a way that, writing u = Uiei+U2e2, the quadratic form Q(u, u) = (Ju\u) 
on E\ assumes the canonical form 

Q(u, u) = nu\ + r 2 u 2 , ri = ±l,r 2 = ±l (1.2) 

Notice that if e\, e 2 is the dual basis, then (jl.2j) means that Jej = Tje*. 
Under these circumstances we want to prove the following 

Theorem 1.1 With the above assumptions and notations, consider the operator 
family H e for e G TZ. Denote: 

H 11 = {H 1 e 1 \ei), H 22 = (#ie 2 |e 2 ), H 12 = (H iei \e 2 ) (1.3) 

Then (ei|ifiei) G 7?., (e 2 |ifie 2 ) G 7?. and taere ezisis e* > such that, for |e| < e*: 

(%) J/n ■ r 2 = — 1, and 

4|# 12 | 2 >(# 11 -# 22 ) 2 (1.4) 
H e has a pair of non real, complex conjugate eigenvalues near A ; 

(ii) If n ■ r 2 = 1 H e has a pair of real eigenvalues near X Q . 
Remarks 

1. The above theorem applies to the PT-symmetric operator family H t = H + 
ieW, where H and iW = Hi above. Here J = P, and hence PH = 
HoP, P(ieW) = -(ieW)P = (ieW)*P so that JH e = H* e J. In that case 
Assumption (jl.4j) follows from the weaker assumption H\ 2 ^ because the 
P— symmetry of Hq and the P— antisymmetry of W entail H\\ = H 22 = 0. 
Indeed, we have Pej = Tjej and 

Hjj = (iW ej \ ej ) = (iPW ej ,P ej ) = {-iWPei\Pe 3 ) = -(iWe^) = -H i5 

2. The physical relevance of Theorem 1.1 is best illustrated by an elementary 
example. Let TC = L 2 (7Z 2 ) and H : 7i — > TC be the (self-adjoint) two dimen- 
sional harmonic oscillator with frequencies u}\ , uo 2 : 

H u = -^An + ^ 2 x 2 + u 2 xl)u 
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We have a(H ) = {E klM } := {kiui + k 2 u 2 + — + —}, h = 0, 1, 2 . . . , % = 1, 2. 
Let again H e = H + ieW, e &7Z, with 

2 

= ^ 9 

v 7 l + xf + xl 

Then is bounded relative to ifo, and PW = —W if Pu(xi, x 2 ) = u(xi, —X2) 
or Pu(xi,x 2 ) = m(— Xi,— x 2 ). Set = l,u 2 = 2, fc x = 2,k 2 = 0; i.e., we 
consider the eigenvalue E 2t0 = E 0j i. Then for |e| > small enough H e has a 
pair of complex conjugate eigenvalues near E 2 $. 

To see this, remark that E 2 $ = E 2 {ui) + E (lu 2 ) = E (ui) + Ei(u 2 ), where 
\i{uji) = (k + l/2)ui are the eigenvalues of the one-dimensional harmonic 
oscillators with frequencies u>i,i = 1,2. E 2 ^ has multiplicity 2. A basis of 
eigenf unctions is given by 

ipi{xi, x 2 ) = e 2 (xi)f Q (x 2 ); ^2(^1, x 2 ) = e (x 1 )f 1 (x 2 ) 

Here eo, e 2 are the eigefunctions corresponding to E (l) and E 2 (l), respec- 
tively; fo, fi are the eigenfunctions corresponding to E {2) and £"1(2), re- 
spectively; note that e , e 2 and f are even while fi is odd. To first order 
perturbation theory, the two eigenvalues Aj(e) : j = 1,2 of H t near E 2t0 are 
given by 

Aj-(e) = E 2>0 + ieXj 
where Xj : j = 1, 2 are the eigenvalues of the 2x2 matrix 

/ Wil^i) \ 

Now is even, ^2 is odd, and W is odd. Therefore T\ ■ r 2 = — 1. Moreover: 
(WV>i|^>i) = (W^lVb) = 0, (W^lV'i) = (^1^2) := w > Therefore 
Aj = ±u> and Aj(e) = £2,0 ±«eu>. Hence the conditions of Theorem 1.1 (i) are 
satisfied and for e small enough H t has a pair complex conjugate eigenvalues 
near E 2t0 . 

3. By essentially the same proof, the result of Theorem 11.11 remains true under 
the following more general conditions: under the above assumptions on H 
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and Hi let Hq admit two real, simple eigenvalues Ex, Ei- Let d :— E2 — Ex be 
their relative distance; D : = dist[(cr(ifo) \ {-^2,-^1}), {-^2,-^1}] their distance 
from the rest of the spectrum; ex, e 2 the corresponding eigenvectors, all other 
notation being the same. Then if d/D is small enough the same conclusion of 
Theorem 11.11 holds provided \eHi 2 \ > 

4. Example: Odd perturbations of quantum mechanical double wells: existence of 
complex eigenvalues. 

d 2 

Let H = L 2 (Jl), H (h) = -h 2 — + x 2 (l + x) 2 , D(H ) = H 2 (K) n L\(K), 

dx l 

W{x) G L% c (ll), \W(x)\ < Ax 4 , \x\ -> 00, W(l-x) = -W{x). Here L\(K) = 
{u G L 2 (1Z) \x 4 u G L 2 (1Z)}. In this case it is known that W is bounded 
relative to H ; moreover d = 0(e' 1/ch ), D = 0(K), w = 0(1) if Ex, E 2 are the 
two lowest eigenvalues, ipx,ip2 the corresponding eigenvectors and w is defined 
as in Point 2 above. Hence the conditions of Theorem 1 are fulfilled in the 
semiclassical regime provided W is continuous at zero with W(0) 7^ and 
that |(ei|We 2 )| > l/C and thus there exist A > 0, B > 0,C > such that 
H e (fi) := H +ieW will have at least a pair of complex conjugate eigenvalues for 

Ae~ B ^ h < ew « Ch. Equivalently, we may consider the double well family 

d 2 2 
Ho(g) = — ~j- 2 +£ 2 (1 + g%) 2 defined on the same domain. Here d = 0(e ' 9 ), 

D = 0(1), w = 0(1). The same argument holds for the general case Ho = 

—h 2 A + V(x), where V : lZ n — ^ 72. is smooth, has two equal quadratic minima 

and diverges positively as \x\ — > 00; W(x) G Lf^ c (JZ n ), \W(x)\ < AV(x) as 

\x\ — > 00 because the estimate for d is the same as above [To"]. 

The second result concerns the opposite situation, a criterion ensuring the reality of 
the spectrum. In this case the natural assumption is the simplicity of the spectrum 
of Hq in addition to its reality. Therefore for the sake of simplicity we assume H 
self- adjoint. 

Theorem 1.2 Let the self-adjoint operator H be bounded below (without loss of 
generality, positive), and let Hx be continuous. Let Hq have discrete spectrum, 
a (Ho) = {0 < \q < \x ■ ■ ■ < \i < ■ ■ ■} , with the property 

5:=w£ [X J+ x-\j]/2>0. (1.5) 
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Assume that all eigenvalues are simple. Then a(H(e)) G 1Z if e 6 TZ, \e\ < . 

II-" ill 

Example 

d 2 

Here again H = L 2 (JZ); H = ~^ + V ( x )^ D ( H o) = H 2 (K)C]D(V). V(x) = kx 2m , 
k > 0, m > 1; W(x) E L°°(K), W(-x) = -W(x). We have: a(H ) = {\ t },n = 
0,1,...; 

A n ~/c 2 ™n m + 1 , n — > oo 

Each eigenvalue A n is simple. Clearly 5 > 1. Denote now if e := Ho + ieiy the 
operator family in L 2 (K) defined by H e = H + H u H l = ieW, D(H e ) = D(H ). 
Then H e has real discrete spectrum for |e| < HM^H^j 1 . 



2 Proof of the results 

Proof of Theorem 1.1 

The proof is based on perturbation theory and consists in two steps. In the first 
one we show that the 2x2 matrix generated by restricting the perturbation H% to 
E\ is antihermitian in case (i) of Theorem 1.1 and Hermitian in case (ii). In the 
second step we show by the method of the Grushin reduction (see, e.g. ^Hl) that for 
e suitably small the control of the above 2x2 matrix is enough to establish the 
result. 

Let {ei,ei} be once more a basis in E\ such that (1.2) holds, and denote by 
e{, e* 2 the dual basis in the dual subspace E* X(i = JE\ . Clearly Jej = Tje*, Tj = ±1. 
We denote n the spectral projection from 7i to E\ . Explicitly: 

n M = (w|ei)ei + (u\e2)e 2 (2.1) 

Consider now the rank 2 operator family U H t Il acting on E\ . The representing 
2x2 matrix is: 

H(e)j, k = X I + eHj ik , Hj >k = (H^e*), j, k = 1, 2 (2.2) 
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Now JH = H*J, m = II* J. We also have JH X = H*J. Therefore: 

(JH^klej) = (Hie k \Jej) = r^H^e*) = TjH) k 
and in the same way 

(JH x e k \e,j) = (HfJe k \ej) = {Je^Hxej) = r^e^i? 1 ^) = T^H^el) = T k E\j 
Summing up: 

Therefore, if T1T2 = 1 the matrix H(e)j >k is hermitian for e G TZ and its eigenvalues 
are real; if instead tit 2 = —1 the matrix H(e)j tk has a real diagonal part and an 
antihermitian off diagonal part for e G TZ and its eigenvalues are complex conjugate. 
This completes the first step. 

We want now to construct an approximate inverse of H e — z near Ao by solving a 
Grushin problem. In this context it is equivalent to the Feshbach reduction, and 
provides a convenient formalism for it. To this end, define the operators 
Vo(z) in the following way: 

R + : H - C 2 , R + u(j) = (u\e*), j = 1, 2; (2.3) 

2 

R- : C 2 -> H, R-U^ = ]T u-(j)ej, (2.4) 
Vo(z) = ( H °- + Z R Q y.VxC 2 ^HxC 2 . (2.5) 

Note that we have identified E Xo with its representative C 2 , and that R + R_ = I, 

the 2x2 identity matrix. 

The associated Grushin system is 

(H - z)u + R^u_ = f 

R+u = f+ 

where u G V, f G H, U-,f + G C 2 . z G C belongs to a neighborhood of A at 
a positive distance from cr(H ) \ {A }. After determining w_ in such a way that 
/ — G (1 — Hq)H the first equation can be solved for u(z) G (1 — n )7Y and 

hence the problem is reduced to the the rank 2 equation R + u(z) = f. To solve 
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explicitly, remark that, for every z in the complex complement of <j(Hq) \ {Ao}, 
Vo(z) has the bounded inverse, 

£n(z)-( E ° {z) E °^ Z) ) (2 7) 

with 



E°_(z) E°_Jz)J : 



E°(z) = (H -z)- 1 (l-U), &4z) = R-, (2.8) 
E°_(z) = R + , E°_ + (z) = (z-\ )I. 

where / is the 2x2 identity matrix. The spectral problem within E\ is thus reduced 
to the inversion of E°_ + (z), and obviously its solution is represented by Ao, eo, e\. 

Now restrict the attention to the set of complex z with dist (z, {Ao}) < 1/2-R, 
where 

R:=\\E°(\ )\\ = \\(l-Tl )(H -\ )- l \\ (2.9) 
so that by the geometrical series expansion 

II^WII * T^x^. < 2 ' 10 > 

Consider the operator from V x C 2 to H defined as 

V e (z) = ( H *- z R ~). (2.11) 



associated to the Grushin system 



(H £ - z)u + R^u^ = f 
R+u = f+ 



(2.12) 



Then 



V e (z)£ (z) = 1 + ^ Q v/ +v/ J=:l + /C. (2.13) 

It is routine to check that P e (z) has the inverse 

£(z)-( Ee{z) E ^ Z) ] (2 14) 



with 



00 ,e 



^) = j2(-) nE °( H i E °) n i ( 2 - 15 ) 



ra=0 Z 
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00 ,e. 



= E(-) n (^i) n£; + ( 2 - 16 ) 

n=0 Z 

oo 

= E(-) ri E°(if 1J E ) n , (2.17) 

n=0 Z 

oo 

= + + Y,(-) nE -( H i E °) n ~ lH i E l- (2-18) 



1 * 

n=l 



where all the series will be proved to have a positive convergence radius (convergence 
means here uniform, or, equivalently, in the norm operator sense). We also recall 
the well known fact that z is an eigenvalue of H € precisely when det Et + (z) = 0. 

We next derive the appropriate symmetries for the inverse operators ^Hl- From 
JH € = H*J we get: 

JR-U- = J2 u -ti) Je i = J2( TU -)ti) e *j> r:= fn 1 r ) 
j=i j=i \ 2 / 

where the second equation follows from 



We thus conclude: 



(R+u\u-) = J2 u -0)( u \ e *j 



JR-U- = R*,tu-, R*_J = tR^ 



Therefore: 



whence 



( J \ ( H e -z R- \ _ ( J{H e - z) JR_ 
\ r J \ R + ) ~ \ tR+ 

(H* - z)J R* + t \ = f (H* - z) R% \ ( J 
R*_J I \ R* J 1 r 



r )W=W(o ° J (2-19) 



Since £(2) = T-'(z) taking right and left inverses we get 



that is 



W J = J W 

EAzY E-JzY Or Or £L(z) E^(z) 



[2.20) 



In particular: 

E^{z)*r = tE_ + (z) 

We can thus conclude that, for z G 1Z, if t± ■ T2 = 1 the 2x2 matrix E-+(z) 
is Hermitian, and antihermitian off the diagonal with real diagonal elements if if 
Ti-t 2 = -1. 

It remains to be proved the norm convergence of the expansions (|2.15|2.17|2.18jl . We 
have, by the relative boundedness condition H-ffi^H < 6||.Ho'0|| + a l|V'll an d Q2.10j) : 

Wh'e'w = H^c^o — - n )II < 

< b\\H Q {H - z)-\i - n )|| + a||(# - - n )|| 

< b\\(H -z)(Ho-z)- 1 (l-U )\\ + 

+ b\z\\\(H - z)-\l - n ) || + a||0Ho - z)-\l - n )|| 

< noil + <k 



for some K{z) > because \z\ < R/2. Therefore 

||£0^l£0)n|| < K n+1^ W^H^ElW < K n+ \ 

WE^H^y'W < K n+ \ \\E°_(H 1 E Q ) n - 1 H l El\\ < K n+1 
Hence the expansions (|2. 1512. 1712. 18)) are norm convergent. 

To conclude the proof we have to verify that the first order truncation of the expan- 
sion for E + (z) yields nonreal eigenvalues, and that the higher order terms can be 
neglected. To this end, first remark that without loss of generality we may assume 
Ao = 0. Then the expansion ()2.18|) yields: 

uniformly with respect to z, \z\ < 1/2R. Therefore 

det£l + (z) = z 2 - (H u + H 22 )tz + e 2 {\H l2 \ 2 + H U H 22 ) + 0(e 3 + e 2 \z\) = 
= [z- e(H n + H 22 )/2} 2 + e 2 [|# 12 | 2 - (H n - H 22 ) 2 /A] + 0(e 3 + e 2 \z\) 

Now detEt + (z), which is real for z G 1Z, clearly has no zeros for z G C, e << \z\ << 
1. On the other hand, for z = 0(e), i.e. z = ew, w = 0(1), 

detE e _ + (z) = e 2 {[w~(Hn + H 22 )/2} 2 + \H 12 \ 2 -(H n ~H 22 ) 2 /4} 
+ 0(e 3 (l + 0(l)) 
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Therefore if 4|ifi2| 2 > (i?n — H22) 2 there cannot be real zeros for e suitably small. 
We can thus conclude that &etE e _ + (z) is zero for z = A±(e), 

A ± (e) = ±[H n + H 22 ± ieyj4\H 12 \*-(H n -H 22 ) 2 ] + 0(e 2 ) 

and this concludes the proof of the Theorem. 



Proof of Theorem 11.21 

Let us first recall that under the present assumptions H t is a type-A holomorphic 
family of operators in the sense of Kato (see Chapter VII. 2) with compact 
resolvents Ve G C. Hence cr(H e ) = {A/(e)} : I — 0, 1, — In particular: 

(i) the eigenvalues A;(e) are locally holomorphic functions of e with only algebraic 
singularities; 

(ii) the eigenvalues A;(e) are stable, namely given any eigenvalue A(e ) of H eo there 
is exactly one eigenvalue A(e) of H t such that lim A(e) = A(e ); 

£—>£0 

(iii) the Rayleigh-Schrodinger perturbation expansion for the eigenprojections and 
the eigenvalues near any eigenvalue A; of H has convergence radius <5//||i/i|| 
where 5i is half the isolation distance of A;. 

Remark that since 5i > 5 V/, all the series will be convergent for all e G Q ro ; 
Q ro := {e G C : |e| < r }, where r := 5/||-£/i|| is a uniform lower bound for all 
convergence radii. 

Assume now without loss of generality, to simplify the notation, ||i?i|| = 1. By 
hypothesis |A; — Az+i| > 25 > 0V7 G Af. First remark that if e G 1Z, |e| < r and 
A(e) is an eigenvalue of H e then |ImA(e)| < 5, i.e. cr(H e ) R C$ = 0, Cs := {z G 
C I |Im z\ > 5}. Set indeed 

R (z) := [H -z]-\ z£a(H ) 

Then \/ z G C such that |Imz| > 5 we have 

\\eHiM*)\\ < M ■ Will ■ \\Mz)\\ < j- n ^77ni < pr^ (2-21) 

dist|z, ct(H )\ \lmz\ 
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Hence the resolvent 



R € {z) := [H e - z}- 1 = R {z)[l + e^Roiz)}- 1 





e 




\lmz\ 


— e 



exists and is bounded if |Imz| > 5 because (I2.21j) entails the uniform norm conver- 
gence of the Neumann expansion for the resolvent: 

oo 

||i? e (^)|| = \\[H e -z)- l \\ = \\Ro{z) £[-effii2o(*)]*|| < 

k=0 

oo 

<l|i?oWIIEl^lll^^)]f< 

fc=0 

Now V/ G Af let Qi(5) denote the open square of side 25 centered at A^. Since 
| A; — Aj + i| > 25, it follows as in ()2.2H) that R e (z) exists and is bounded for z G dQi{5), 
the boundary of Qi(5). We can therefore, according to the standard procedure (see 
e.g.jTT], Chapter III. 2) define the strong Riemann integrals 

Pi(t) = 7T^f Re(z)dz, 1 = 1,2, 
2m JdQi{&) 



As is well known, Pi is the spectral projection onto the part of cr(H e ) inside Q[. 
Since H e is a holomorphic family in e, by well known results (see e.g. Thm. 
VII.2.1), the same is true for Pi(e) for all I G A/". In particular this entails the 
continuity of Pj(e) for |e| < r . Now Pi(0) is a one-dimensional: hence the same is 
true for Pi(e). As a consequence, there is one and only one point of cr(H e ) inside 
any Qi. Now cr(H e ) is discrete, and thus any such point is an eigenvalue; moreover, 
any such point is real for e real because <j(H e ) is symmetric with respect to the real 

oo 

axis. Finally, we note that if z G 71, z ^ [J] A; — 5, A/ + 5[ the Neumann series (I2.21J) 

i=i 

is convergent and the resolvent R e (z) is there continuous. This concludes the proof 
of Theorem 1.2. 
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